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Abstract:  This  paper  adresses  the  problem  of  assesing  the  quality  of  an  a-posteriori 
error  estimate  of  a  finite  element  solution.  An  error  estimate  based  on  local  L2-projections 
is  analyzed  in  the  case  of  translation  invariant  meshes.  It  is  shown  that  for  general  meshes 
this  technique  does  not  lead  to  an  asymptotically  exact  estimator.  The  problem  is  analyzed 
in  detail  in  the  one-dimensional  setting.  It  is  shown  that  an  asymptotically  exact  estimator 
is  not  the  optimal  one  when  the  solution  is  not  sufficiently  smooth.  An  optimal  estimator 
for  adaptively  constructed  meshes  is  given.  Finally,  a  general  mathematical  framework  for 
the  quality  assesment  of  estimators  is  introduced. 
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1.  Introduction.  Since  the  first  papers  by  Babu§ka  and  Rheinboldt  [5,6]  on  a-poste- 
riori  error  estimates  in  the  finite  element  method,  the  subject  has  become  increasingly 
important  in  finite  element  practice  [9,10,12,14,16,19,20,23,25,30,32].  Presently  there  exist 
various  codes  (for  research  or  for  commercial  use)  including  a-posteriori  error  estimation. 

There  are  essentially  three  major  types  of  estimators  (see  also  the  discussion  in  [27]): 

a)  estimators  based  on  residual  considerations, 

b)  estimators  based  on  averaging  techniques, 

c)  estimators  based  on  extrapolation. 

The  first  proposed  estimators  were  of  residual  type.  For  some  of  them  a  rather  precise 
mathematical  analysis  has  been  presented  [1,4,11,17,18,23].  Others  were  derived  on  purely 
heuristic  grounds.  Estimators  of  this  type  have  been  also  used  as  the  basis  for  adaptive 
approaches  (see  e.g.  [8,16,22,23,25]).  The  first  research  adaptive  code  based  on  a  residual 
estimator  is  likely  FEARS  [7,20]. 

In  recent  years,  the  estimators  of  type  (b)  have  become  extremely  popular;  in  particu¬ 
lar,  the  estimator  proposed  by  Zienkiewicz  and  Zhu  [31].  This  paper  will  address  estimators 
of  this  type  in  detail. 

Estimators  of  type  (c)  are  usually  used  in  connection  with  the  p-version  of  FEM  (see 
e.g.  [26,27]),  while  estimators  of  type  (a)  and  (b)  are  used  mostly  but  not  exclusively  in 
connection  with  the  h-version  of  the  method. 

2.  The  model  problems.  We  shall  consider  the  two  following  model  problems. 

a)  Let  ft  C  R2  be  a  bounded  polygonal  domain  with  boundary  T.  We  assume  that 
r  =  Td  U  r.  (for  simplicity  we  assume  that  Td  has  positive  length).  We  are  interested  in 
the  solution  of  the  problem 

{-Au  =  /  ,  in  ft  , 
t»  =  0  ,  on  I'd  , 
fi=?,  on  T.  ; 

we  denote  by  n  the  outer  normal  vector  to  I\  We  assume  that  /  and  g  are  sufficiently 
smooth  functions;  then  the  solution  of  (2.1)  exists  and  is  unique. 

b)  Let  ft,  :=  (0,1)  x  (0,1);  we  consider  the  Laplace  equation  with  periodic  boundary 
conditions: 

-Au  —  f  ,  in  ft,  , 

«(0,x2)  =  u(l,x2)  ,  ^(0,x2)  =  $j-(1,x2)  ,  0  <  x2  <  1  , 

u(xi,0)  =s  u(xi,l)  ,  ^-(xj.O)  =  ^(xi,  1)  ,  0  <  xi  <  1  , 

Jn  u  dx\dxi  =  0  . 

We  assume  that  /  is  sufficiently  smooth  and  that  /n<  /  dx\dx?  =  0;  then  the  solution  of 
(2.2)  exists  and  is  unique.  If  /  and  u  are  periodically  extended  to  R2  then  obviously  u  solves 


3 


the  differ'ntial  equation  in  RJ  also.  From  now  on,  when  we  speak  about  the  regularity  of 
the  solution  of  this  problem,  we  mean  the  regularity  of  the  periodically  extended  solution. 


Remark  2.1.  We  restrict  ourselves  to  these  model  problems  for  the  sake  of  simplicity 
but  our  results  hold  in  more  general  situations.  □ 

Let  |u|j  :=  /n  |Vu|2diidi2  and  |u|2  :=  fn  |u|2dxjdx2;  let  :=*  {u  :  |u|2  -t  |u|2  < 

00}  be  the  usual  Sobolev  space.  Further,  let  us  define  ifg(Q)  :=  {u  £  H'(Q)  u  = 
0  ,  on  Td}  and  #JER(^>)  :=  {u  €  u  has  period  1  and  /n<  u  dx\dxi  =  0}. 

Obviously  |  •  |j  is  a  norm  on  Hq(£1)  and  on  H*KR(Q,). 

We  consider  the  usual  variational  form  of  (2.1)  and  (2.2).  Let 

(2.3)  2?(u,t;):=  f  Vu  •  Vv  dxidx2 

Jo 

and 

(2.4)  H|2  :=B(«,u)*|u|J  . 

Further,  for  problem  (2.1),  let 

i2.5)  F(v)  :=  I  fv  dx\di2  +  f  gv  da 

Jo  J  r„ 

and  for  problem  (2.2)  let 

(2.6)  F(v)  :=  f  fv  dx\dxi  ; 

Jom 

both  are  linear  functionals  on  Hq(Q)  and  respectively.  The  weak  solution  uo  € 

(respectively  «o  €  #RER(fti))  is  such  that 

(2.7)  B(u0,v)  =  F(v)  ,  VveffdW  (resp .  v  €  ^(R(fi,))  . 


For  the  finite  element  method  we  partition  Q  into  a  set  of  (dosed)  elements  r  defined 
by  the  mesh  7*;  (to  fix  ideas,  let  us  consider  a  triangular  mesh).  Let  the  mesh  be  regular 
in  the  sense  of  a  minimal  angle  condition  and  let  h  denote  the  maxr gr*  diam(r). 
Furthermore  let 


S*r(fi)  :=  {u  £  Hl(Q)  :  u|r  €  Vr  ,  Vr  €  Th)  , 

{Vr  denote  the  set  of  polynomials  of  degree  r)  and  let 

s;,0(ft)  :=  SI(ft)  n  ff‘(fi)  and  5Ar,p„(ft,)  :=  Sr,(Q)  n  H'rtR(Q,)  . 
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The  functions  of  S£FER(fl,)  can  be  periodically  extended  to  R2;  we  denote  SJ  pER  to  the 
set  of  these  extensions. 

We  say  that  a  mesh  is  (hi, /^-translation  invariant  (hi,hj  positive  numbers)  if  for 
any  pair  of  integers  V\,V2  and  for  any  u  €  SJ>pER,  u(xi  +  i'iAi,x2  +  1^2)  6  SJ  pER.  Figure 
2.1  shows  some  of  these  meshes. 


(b) 


Figure  2.1 
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a)  (h,  A)- translation  invariant  mesh, 

b)  (2h,  2A)-translation  invariant  mesh, 

c)  (2 h,  A)- translation  invariant  mesh. 

The  finite  element  solution  uFE  €  SJ  0(il)  (resp.  SJ^pER(fl,))  is  defined  as: 


(2.8)  B(un,v)  =  F(»)  ,  Vu  €  SJ(0(fi)  (resp.  v  €  Sjk>PER(fl,))  • 


3.  A-posteriori  error  estimates  based  on  averaging  techniques.  The  main 
idea  behind  the  averaging  techniques  is  to  observe  that  VuFE  £  5£(fi)2  and  to  construct 
U  €  SJ(fl)2  using  only  urE  with  the  hope  that 


(3.1)  |Vu0 -U|0  <  |Vu0  -  VuFe|0  . 
If  (3.1)  holds  then 

(3.2)  |U  -  VuFE  j0  ss  ||u0  -  uFE  ||,  • 
and  we  can  define 


(3.3)  e  :=  |U  -  VuFE|0 

as  the  a-posteriori  error  estimator  (with  respect  to  the  energy  norm). 

Given  a  class  of  functions  H  C  Bq(Q)  (respectively  H  C  /fpER(fl, ))  we  call  £  a  correct 
estimator  on  H  if  there  exist  constants  C\  and  C2  (0  <  C\  <  C2  <  00)  such  that 


(3.4) 


Cx  < 


u0  -  u 


<C2  ; 


r*lls 
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whenever  the  solution  uq  belongs  to  H\  these  constants  must  be  independent  of  the  mesh- 
size  but  could  depend  on  the  minimal  angle  of  the  mesh  and  on  the  class  H. 

On  the  other  hand,  the  estimator  e  is  called  asymptotically  exact  on  H  if 


lluo  - 


1  ,  as  h  — *  0  , 


whenever  the  solution  uo  belongs  to  H. 

There  are  many  suggestions  in  the  literature  for  the  construction  of  U  (and  hence  of 
e),  mostly  without  any  theoretical  justification.  Let  us  divide  these  constructions  into  a 
few  basic  groups: 

1)  global  construction  of  U, 

2)  local  construction  of  U, 

3)  semilocal  construction  of  U. 

Further  we  may  distinguish  between  two  other  groups: 

a)  general  construction  of  U, 

b)  construction  of  U  based  on  the  available  information  of  the  problem  (i.e.  the 
differential  operator)  under  consideration. 

In  what  follows  we  restrict  ourselves  mostly  to  the  case  r  =  1  (i.e.  piecewise  linear 
elements)  for  simplicity.  In  this  case  the  solution  uFE  and  U  itself  are  characterized  by 
their  nodal  values.  Each  nodal  value  of  U  depends  only  on  the  values  of  uFE  in  some 
neighborhood  of  this  nodal  point.  Different  kinds  of  neighborhoods  distinguish  the  groups 
1,  2  and  3  mentioned  above. 

Let  A  be  a  nodal  point  of  a  mesh  7*  and  s  >  1  an  integer.  For  s  =  1  we  define 

Afl(^):=U{r€T*  :  r  9  A} 

and  for  s  >  1  we  define  recursively 

W"+1(A)  :=  1J  M\B)  . 

B  node  of  Af’(A) 

When  U(A)  is  defined  only  by  the  values  of  uFE  on  we  say  that  the  construction 

is  local.  When  the  values  on  N*(A)  for  s  >  1  are  used,  we  say  that  it  is  semilocal.  The 
construction  is  global  if  the  values  of  uFE  on  the  whole  domain  Q  are  used.  Figure  3.1 
illustrates  these  notions. 
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o  nodal  points  of  M"1 

Figure  3.1 


A  typical  construction  of  U  :=  {U\,Ui)  consists  of  using  the  L2(ft)-projection  of  • 
on  SJ(fi)  as  Ui,  i  =  1,2  (see  e.g.  [31]).  This  construction  belongs  to  the  category  1-a;  the 
matrix  for  determining  the  nodal  values  of  Ui  is  the  usual  mass  matrix.  A  local  construction 
of  type  2-a  is  obtained  if  the  lumped  mass  matrix  is  used  instead  of  the  consistent  one. 


U,(A)  can  also  be  defined  as  the  value  at  that  node  of  the  L2  ( AT1  ( .4)  )-pro jection  of 
onto  the  set  of  continuous  piecewise  linear  functions  on  that  patch.  This  construction 
also  belongs  to  category  2-a  and  will  be  especially  addressed  later. 


In  general  any  local  construction  defines  U,{A)  as  a  weighted  average  of  the  constant 
values  of  -gff-  in  the  elements  t  C  Af'iA).  For  various  types  of  averages  we  refer  to  [23,29]. 

These  kind  of  constructions  usually  lead  to  estimators  correct  on  Hq(Q)  (respectively 
on  ffpER(^»))  in  the  sense  of  (3.4)  with  constants  C\,Ci  depending  on  the  geometry  of 
the  mesh.  However,  asymptotic  exactness  is  related  with  superconvergence  effects.  If  the 
mesh  is  translation  invariant  then  many  superconvergence  results  hold;  we  refer  to  section 
30  of  [28]  for  this  subject.  These  results  can  be  easily  employed  for  the  construction  of  U 
and  hence  of  the  estimator. 


Let  us  now  consider  problem  (2.2).  Assume  that  the  mesh  is  (h\,  /^-translation 
invariant.  Let 


(£},*«)  (*1>*2)  :=  jT-  u(xj  +  y,I2)  -  u(r,  -  y,I2)  , 

(D\,ku)(x u(l|,X2  +  y)  -  u(li,X2  -  y)  . 

Then  we  have  (see  Theorem  30.1  of  [28]): 
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THEOREM  3.1.  If  the  periodic  extension  of  the  solution  uo  of  (2.2)  has  three  bounded 
derivatives,  i.e.: 


(3.6) 


d3 


uo 


-(Xi ,  x2) 


|  dxidxjdxk 

and  the  mesh  is  (hi, h2)-translation  invariant,  then 


<M,  V(xi,x2)  €  R2  ,  i,j,k=  1,2, 


( 

\dXi 


-  Dhu 


F* 


<  C\\ogh\h2M  , 


*-1,2. 


Remark  3.1.  Theorem  3.1  can  be  easily  understood.  is  the  solution  of  problem 
(2.2)  with  /  replaced  by  and  D]k uFe  is  the  finite  element  solution  of  problem  (2.2) 

with  /  replaced  by  D\hf.  Since  —  f?,1  h/j  =  0(h2),  then  the  conclusion  of  the  theorem 
holds.  (The  log(/i)-term  comes  from  the  L°°  error  estimate  of  the  finite  element  solution). 
□ 

This  theorem  allows  us  to  design  proper  averages,  for  n  s  such  as  those  in  Fig.  2.1. 
For  the  following  examples  let  us  assume  that  the  solution  uo  of  problem  (2.2)  is  in  the 
subspace  H  of  periodic  functions  with  three  bounded  derivatives. 

Mesh  2.1. a.  Consider  the  set  AT1  (A)  shown  in  Fig  3.2. 


—  h  — h  — 


Figure  3.2 


We  have: 


+  \'og h\o(h2) , 
~(M2)  =  ~f(M2)  +  |logfi|<9(h2)  , 
|^W)  =  |^(ivi)  +  |IoS^lC?(fi2)  , 
| ^W)  =  ~^-(N2)  +  \\ogh\0(h2)  . 
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The  average  of  on  the  right  hand  side  of  (3.7)  is  identical  to  the  nodal  value  Ui{A) 
obtained  by  the  I2(iV1(.4))-projection  (or  equivalently  by  the  lumped  £2(fi)-projection). 
An  analogous  result  is  true  for  U?.  Therefore,  since  U  is  defined  by  its  nodal  values  -.'  hich 
are  higher  order  approximations  of  the  exact  solution,  an  explicit  calculation  shows  that 
(3.1)  holds  and  hence,  on  these  kind  of  meshes,  this  projection  leads  to  an  asymptotically 
exact  estimator  on  the  space  H  of  periodic  functions  satisfying  (3.6). 

Mesh  2.1. b.  In  this  case  the  mesh  is  (2h,  2/i)-translation  invariant.  Consider  A’1  (-4) 
and  Af'{C)  as  shown  in  Fig  3.3. 


£>4  JV,  Z?3  Ns  £>2 

T 

h 

Br  |  Ds 

h 

i 

De  Ni  D7  N4  D9 
b  )AT(C) 

Figure  3.3 

We  have 

&U.  i 

g^(A)  =  ^  -  U-(^3)]  +  \\ogh\0(h2) 

=  ±  {[«„(*,)  -  u„(A)J  +  (upb(A)  -  ufe(53)]}  +  | log/i|C>(/i2) 

(3S)  =\  %7{M2)  +  lh7{M,)  -Hlo«W*2)- 


Bi 


B< 

*—  h  — ►«—  h  —* 

a)  AH  (A) 
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Similarly  for  A  1  (C): 


dz,( ]  ~  8 


ax,  <'v')  +  2tf('V!)  +  %f('v>,  +  lff!'v') 


(3.9) 


du 


■^(NS)  + 

ox  i 


+  i  log  *|<P(A2) 


It  is  easy  to  check  that  the  L2(A’1)- projection  of  gives  the  same  averages  as  in 
(3.8)  and  (3.9).  Therefore,  also  for  these  kind  of  meshes,  this  projection  leads  to  asymptot¬ 
ically  exact  estimators  on  the  space  H  of  periodic  function  with  three  bounded  derivatives. 

Mesh  2.1. c.  Now  the  mesh  is  (2h, /^-translation  invariant.  Let  .'/’(A)  be  as  in  Fig  3.4. 


Figure  3.4 


In  this  case,  the  value  obtained  by  the  £2(A’1  )-projection  is 


(3.10) 


<MA)  =  - 


1 2%7^+2iZw+%7^+d-£r^\ 


This  value  coincides  up  to  higher  order  terms  with  fj^(C)  (see  Fig.  3.3;  dist(AC)  =  y) 
instead  of  f^-(A).  An  analogous  statement  is  true  for  each  interior  node;  for  instance. 
Ui{Bi)  coincides  with  § ^{D)  (see  again  Fig.  3.3;  dist(Si£>)  =  y)  instead  of 

Let  R  be  the  square  of  vertices  A,  B\ ,  P,  B^\  by  usmg  the  corresponding  averages  of  the 
form  (3.10)  for  each  vertex  and  applying  Theorem  3.1  and  appropiate  mean  value  theorems 
we  obtain: 


L 


dup 

dxi 


2 

diidi2 


h*_ 

27 


d2up 

dx\dxi 


2 


4-  |logh|0(h5)  . 


Therefore,  U\  is  not  an  approximation  of  of  higher  order  than  A  similar  analysis 

is  valid  for  Up..  Hence  the  condition  (3.1)  on  which  the  effectiveness  of  this  estimator  is 
based  is  not  satisfied. 
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So,  the  estimator  t  will  not  be  asymptotically  exact  on  these  kind  of  meshes  unless 
a  fortunate  compensation  happens.  But,  in  general,  this  is  not  the  case.  For  instance 
consider  the  following  problem: 

f  -  Au  as  0  ,  in  fi  , 

1  u  =  9  ,  on  T  , 

where  fl  is  the  unit  square  and  g  is  defined  in  Figure  3.5: 


*2 


The  solution  of  this  problem  is  uo(  11,12)  =  X1I2.  In  this  example  the  F.E.  approxima¬ 
tion  coincides  with  the  exact  solution  at  the  nodes  of  the  mesh  (i.e.  uFe  is  the  interpolant 
of  uq  ).  So  the  exact  error  and  the  estimator  can  be  explicitly  computed.  For  any  triangle 
r  of  the  mesh  we  have 


/ 


2  h 4  r 

|Vu0  -  VuFC|  =  —  and  j 


~/i4 


iU  -  V“'*'  *  IT 


Therefore 


i  :=  r. - - - —  =  ~  *  0.88  , 


independently  of  the  mesh  size  h.  That  is,  in  spite  of  the  smoothness  of  the  solution,  £  /♦  1 
as  h  — »  0. 

This  example  shows  that  the  estimator  based  on  the  L2(.V]  )-projection  is  not  always 
asymptotically  exact,  even  in  the  simple  case  of  translation  invariant  meshes.  To  obtain 
an  asymptotically  exact  estimator,  the  geometry  of  the  mesh  should  be  used  in  a  more 
sophisticated  way. 

Now,  let  us  briefly  consider  the  case  r  =  2.  For  (hi ,  /^-invariant  translation  meshes, 
let 

For  smooth  functions  u  this  is  a  fourth  order  approximation  of  The  following  theorem 
is  also  a  consequence  of  Theorem  30.1  of  [28]. 
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THEOREM  3.2.  If  the  periodic  extension  of  the  solution  uq  of  (2.2)  has  four  bounded 
derivatives,  i.e.: 


d*uo 


■(xux2) 


<  M  ,  V(x1,x2)€Rj,  i,j,k,l  =1,2, 


|  dxidxjdxkdn 
and  the  mesh  is  (hi ,  h2)~ translation  invariant,  then 

<ChHf,  t  =  l,2. 


(ax!  DUu**)  (*''**) 


By  using  this  theorem  we  may  construct  U  €  S\(Q)2  satisfying  (3.1)  for  invariant 
translation  meshes.  For  instance  consider  a  uniform  mesh  as  that  in  2.1. a. 


h 


Figure  3.6 


We  need  to  define  the  nodal  values  of  U.  For  instance,  with  the  notation  of  Fig.  3.6, 
let 

U\(C)  :=  (D\tkun)(C) 

4ure(-P)  ~  ur*(B)  _  1  urz(E)  ~  “fJ-4) 

3  h  3  2  h 

2 


r**re(0)  -  »f*(Q  +  -  ufe(B) 


hi  2 


h/2 
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\urt.(E)  -  un(D)  «„(£)  -  uri(C)  ,  urz(C)  ~  urz(B)  -  u„(A) 

•  I  '  I  ' 


h/2 


by  using  Theorem  3.2  w»  have 


A/2  A/2 

^  f^-(P)  +  ^(A/) 

12  ax1  axj 


A/2 


^(C)  =  |^(C)  +  CJ(A3). 

The  same  construction  can  be  made  for  any  other  node  P  of  the  mesh  and  also  for  U2\ 
i.e.,  in  general, 

U,(P,:=(DlhurE)(P). 
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It  can  be  easily  proved  that  this  construction  yields  an  asymptotically  exact  estimator. 

We  observe  that  in  this  case  the  construction  of  U  is  semilocal;  in  fact  U(C)  is  an 
average  of  -jjf*-,  i  —  1,2  at  points  of  triangles  r  C 

In  the  case  of  problem  (2.1)  we  have  to  extend  the  solutions  uo  and  uFE  outside  of  ft 
so  that  the  averages  can  be  computed  also  on  the  boundary.  We  suggest  the  use  of  Babic 
extension  techniques.  In  the  next  section  we  shall  discuss  in  detail  these  techniques  for  the 
one-dimens;onal  case.  For  the  n-dimensional  case  see,  for  instance,  [21]  (pp  75-77). 

For  general  meshes  it  is  probably  impossible  that  estimators  of  this  kind  can  be  asymp¬ 
totically  exact,  since  this  seems  to  depend  on  a  superconvergence  effect.  At  any  rate,  for 
practical  constructions,  those  estimators  that  are  asymptotically  exact  for  invariant  trans¬ 
lation  meshes  seems  to  be  preferable. 

We  addressed  only  local  and  semilocal  constructions.  However,  it  is  very  likely  that 
global  estimators  based  on  X2-projections  will  not  bring  a  larger  benefit  and  will  not  avoid 
the  problem  with  general  meshes. 

So  far  we  have  discussed  only  estimators  of  type  (a).  Semilocal  or  globad  constructions 
of  U  of  type  (b)  are  more  expensive  but  could  lead  to  asymptotically  exact  estimators  also 
for  general  meshes.  These  constructions  are  based  on  postprocessing  techniques,  but  we 
will  not  analyze  these  kind  of  estimators  here.  For  further  references  on  this  subject  see 
[15,23], 

Remark  3.2.  It  is  very  likely  that  the  first  result  related  to  the  superconvergence  of 
smoothening  constructions  that  could  be  used  to  define  a-posteriori  estimators  is  in  [24] 
(1969).  An  important  result  of  the  same  type  is  also  in  [13].  □ 


4.  Model  problem  in  one  dimension.  In  the  forthcoming  sections  we  shall  analyze 
in  detail  the  one  dimensional  case  to  gain  additional  insight.  Let  us  consider  the  problem 

f  -u"  =  /  ,  on  /  :=  (0, 1)  , 

(41)  {  u(0)  =  u(l)  s  0 


and  let  /  be  such  that  u  €  H<j(/). 

Let  T  be  a  partition  of  the  interval  I  with  nodes  0  =  Jo  <  xj  <  •••  <  x„  =  1; 
:=  [xi_i ,  Z|],  h,  :=  x,  —  x,_j ,  i  —  1, . . . ,  n. 

Consider  the  finite  element  method  with  piecewise  linear  elements  on  the  mesh  T.  In 
this  case  we  have 


uPE(x,)  =  “(*.)  >  »  =  1 . n  : 

i.e.  uFB  is  the  piecewise  linear  interpolant  of  u. 

Denote  t,  :=  u  -  uFB.  We  shall  be  interested  in  error  indicators  q,  approximating 


Klu. :»  KUj,  <-  [/,,«)’* 

(4.2) 


,  where 


1?  :=  /  (tf-O2; 
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here  U  is  a  continuous  piecewise  linear  function  on  I.  In  Fig.  4.1  we  show  u',  u'Ft  and 

U\i,. 


We  shall  consider  only  local  constructions  of  U  and  discuss  various  aspects  of  such 
constructions.  For  any  x„  1  <  t  <  n  —  1,  we  define 

>-t+i 

(4-3)  U(xi)  :=  Y,  7j°ttw(*j)  . 

j*i-i 


where  7*,)  depend  on  the  mesh. 

We  constrain  the  coefficients  7*1*  so  that  U{xt)  =  u'FI(x>)  when  u'FE  is  linear  on 
^(i,):=/.U/,+1.  This  leads  to  the  formula 


U(zi)  = 


hi 


+  a, 


^Fe(^»)  UFE(Xj_l) 


‘i+1 


hi 


(4.4) 


=  “p«k  +  KJ,  -  «'„!/!♦.  +  Pi  KJ,  , 


where 


(4-5)  K«],  :=«p,I/<+.  -“'„k 

is  the  jump  of  the  derivative  of  at  n\  (obviously  a<  =  1  +  Pi).  U(xi)  depends  only  on 
the  values  of  u'FE  in  the  elements  I,  and  Ii+i  which  contains  the  mesh  point  x,. 

For  the  values  of  U{xo)  and  U{xn)  formula  (4.4)  cannot  be  directly  applied.  To  use  this 
formula  we  need  to  extend  u  and  uFE  for  x  <  0  (and  for  x  >  1)  by  a  process  that  preserves 
the  smoothness  of  the  function  under  consideration.  For  example,  it  is  possible  to  use 
BabiC  extension  (see  e.g.  (21),  pp  75-77).  This  extension  preserves  the  smoothness  of  u  in 
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the  sense  that  if  u  €  ffJ((0,ij]),  (1  <  s  <  3),  then  the  extended  function  u  on  [x_i,X2] 
which  coincides  with  u  on  (xo,xjJ  is  such  that  v  6  /f,([x_i,X2])  and  IMI/f  *({*_,,*,])  < 
C||u||w,([Xo  Zjj),  where  C  depends  only  on  the  local  regularity  of  the  mesh.  Using  Xi  *  -/ij, 
u(x_j )  depends  only  on  u(xo),  u(xi )  and  u(x-N  and  hence  the  extension  can  be  used  for 
uFB  as  well.  When  1  <  s  <  2  the  simple  antisymmetric  extension  (that  is  a  BabiC  extension 
preserving  at  most  H2- smoothness)  can  be  ur  :i.  In  any  case,  (4.4)  can  be  applied  to  the 
extended  function.  Let  us  remark  that  if  u  is  e:  tended  antisymmetrically  U(x0)  =  u'FE(x0 ). 

Now,  defining 

(4.6)  V(x):=U(x)-u'n(x),  x€/, 


we  have  for 


(4.7)  V(x)  =  V,(x)  :=  0i.i  Vi-,(x)  +  a,  KJ,  *,(x)  , 

where  and  ^i-i  are  the  standard  linear  shape  functions  on  /,.  Hence  using  (4.2)  we  get 

(4.8)  tf-MI Ji[Vx{x)\2dx 


and  so  we  can 

KJ,;ie-: 

(4.9) 


Vi  =  *  J  \vi\xn  ax 

understand  the  indicator  rjf  as  a  quadratic  form  on  the  jumps  [u'FE]  ,_j  and 


?  =  k>  £  u*Ke]>  K«l*  • 


Using  (4.7)  and  (4.8)  7,,*  are  directly  related  to  &•_  j  and  or,  in  an  obvious  way.  However 
the  indicator  (4.9)  can  be  defined  without  any  relation  to  the  function  VJ(x).  We  can  simply 
define  and  address  the  accuracy  of  that  indicator. 

We  define  the  error  estimator: 

(4.io) 

i=i 


Let  us  mention  that  the  well  known  Zienkiewicz-Zhu  estimator  (Z-Z)  [31]  which  uses 
the  L2-projection  of  u'FE  on  ;V*(x,)  leads  to 

/,  ,,,  a  *»«-i  .  ^*+i 

(4.11)  0i-l  —  — T  1  >  a<  =  L  ,  L 

« t-1  +  hi  "i  +  «i+l 

On  the  other  hand,  the  Babuska-Miller  estimator  [4]  (B-M)  leads  to 

(4'12)  /J”‘  =  “hi^+hi  ’  01  =  hi  +  A,+1  ‘ 

Obviously  for  uniform  meshes  both  estimators,  Z-Z  and  B-M,  are  identical. 


(4.12) 


A- 1  =  - 


A-i  =  - 


By  using  ficticious  nodes  xj  :=  — /»2  and  xn+1  :=  1  +  /»n_i,  the  Babic  extension  of  u 
preserving  H‘ -smoothness  for  1  <  s  <  3,  gives  [u'Ft]0  =  [u'FE],  and  [u'FJn  =  [u'FE]n_,. 
The  Z-Z  estimator  uses  [«FB]0  =  (uFl]n  =  0  (these  are  also  the  values  obtained  by  the 
antisymmetric  Babic  extension).  Nevertheless  it  is  easy  to  see  that  both  extensions  lead 
to  the  same  indicator. 
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5.  Principles  of  the  assesment  of  the  quality  of  an  error  estimator.  The 
formulation  of  the  quality  of  an  error  estimator  is  far  from  being  obvious.  One  of  the  most 
natural  principles  is  to  compare  the  indicator  and  the  true  error  |eM |j  To  this  end  we 
define  the  elemental  effectivity  indices 

(5.1)  &  :=  —  ,  i  =  1, . . . ,  n  . 

We  say  that  the  indicator  is  asymptotically  exact  if  for  any  mesh  such  that  0  <  C\  < 
X7^7  <  Ci  <  oo,  (x  *  1, . . .  ,n  —  1),  and,  for  u  sufficiently  smooth  and  such  that  |e,|j  {  > 

Chj ,  ( C  >  0),  we  have  for  all  x  =  1, . . . ,  n, 

Zi  -*  1  ,  as  hi-*  0  . 


If  these  elemental  effectivity  indices  are  used  as  a  criterion  to  choose  an  indicator,  the 
asymptotically  exact  indicators  would  be  preferable  when  the  mesh  is  sufficiently  refined. 
However  we  note  that  this  obvious  criterion  needs  to  be  refined  as  it  will  be  seen  later. 

In  practice,  the  solution  u  is  unknown  but,  usually,  it  is  known  that  u  belongs  to  certain 
space  of  functions  7i  C  Hq  (/).  In  this  case  we  may  define 


(5.2)  (£«)(*)  :=  (C,u)( x)  :=  (e'.  -  V,)(x)  ,  x  €  /.  , 

and  understand  C,  as  a  mapping  of  H  onto  the  space  L2(It).  Since 


(5.3) 


W  ~  ~  IKIo./.  -  le*  ~  V.lo,A  —  > 


then  |£,u|0  7j  can  be  used  for  the  assesment  of  the  quality  of  the  indicator. 

We  could  introduce 

p{7i)  :=  sup  -r— - 

as  a  measure  of  this  quality.  For  any  indicator,  there  exists  a  function  v  €  H  such  that 


|£.«H0,/,  =  P(^)|«*  |if/. 


If  we  do  not  know  anything  more  about  the  solution  than  ugW,  the  safest  choice  would 
be  to  use  that  indicator  which  minimizes  p^H).  However  this  approach  has  a  major  disad¬ 
vantage:  f^H)  <  oo  only  for  very  restrictive  spaces  H.  For  example,  even  for  H  =  C°°(/), 
p(H)  =  oo 

To  avoid  this,  we  can  asses  the  quality  of  an  indicator  by  the  norm  of  C, 


(5.4) 


IIAII  -•=  sup  ^  J07, 

INI* 
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and  prefer  those  indicators  leading  to  smaller  j|£,||.  This  quality  measure  obviously  de¬ 
pends  on  the  space  H(7)  of  possible  solutions.  In  subsequent  sections  we  shall  study  this 
criterion  in  detail. 

So  far,  we  described  criterions  of  elemental  type.  We  can  also  define  the  quality  of  the 
indicator  through  the  corresponding  estimator,  i.e.  to  use  the  global  effectivity  index 


(5.5) 


gr-??)* 


or  analogously  the  norm  of  £ 


(5.6) 


]£<,  [//(e',-r)a 
11  11  ‘  I'Si  INI*  '  IMI* 


and  to  prefer  an  indicator  based  on  the  (global)  quality  assesment  of  the  estimator. 

6.  Comparison  of  Z-Z  and  B-M  indicators  based  on  the  elemental  effectivity 
index.  In  this  section  we  compare  Z-Z  and  B-M  indicators  in  order  to  show  the  difficulties 
arising  in  the  comparison  of  any  estimators. 

THEOREM  6.1.  B-M  indicator  is  asymptotically  exact. 

Proof.  It  can  be  easily  seen  that  when  u  is  a  quadratic  polynomial,  U  coincides  exactly 
with  u'  and  hence  =  1.  In  general,  for  n  sufficiently  smooth  (say,  for  instance,  that  u"'(x) 
is  be<mded  for  x  €  /<,  where  I,  :=  /,_i  U  /,  U  /,+  i),  let 


u(x)  =  Co  +  Cj(l  -  Xj)  +  c2(x  -  x,)2  +  R{x) 


be  the  Taylor  expansion  of  u  and  let 

u(x)  :=  Co  +  ci(x  -  x,)  +  c2(z  -  x,)2  . 

Let  u„  be  the  finite  element  aproximation  of  U  (i.e.  its  interpolant).  Let  U  and  U  be  the 
linear  functions  on  /<  defined  by  (4.4)  by  using  B-M  coefficients  (4.12)  for  uPE  and  uFE 
respectively;  i.e.,  for  j  =  i  -  1,*, 


U{*i)  = 

a(xi+i) -u(x>)  , 

hj+i 

u{x})  -  u(lj.i) 

h, 

+  ^>+1 

hj  +  hj+i 

h, 

U{x})  - 

u(xi+1)-u(ij) 

'Vj+l 

«(*>)-  “(xj-i ) 

*7 

+  ^»;+l 

hj+i 

hj  +  hj+i 

h, 
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Since  u  is  a  quadratic  polynomial,  then  U(x)  *  u'(x),  for  x  €  I  and  hence 


*  -  l«.U|  =  ||P  -  - 1“'  -  |  SIP- «'Im,  SIP- PI,.,.  ♦I*'- »%./, ; 

therefore. 


(6.2) 


16  -11  = 


Vi 


-1 


^  lU-Ulojt+W-u'lo 


Now,  since  u  —  u  =  f?,  then  by  using  (6.1)  we  have  for  j  *  i  —  1,  i, 

in  frx,_  x  f?(xJ+1)  .  h;+i  fl(x,)  -  fl(i>_i) 

+  ft-  +  *<+1  h, 


_ a>t  /■  \  . 


+  fy+l 


#(C;+l)  + 


ff(Ci) 


with  Ij_!  <  C>  <  Xj  <  Cj+1  <  Xj+i;  hence 

|tf-tf|0/i<h*max|#(x)|. 

<€/i 


On  the  other  hand, 


Iu'-^Va  1  niax|f?'(r)| 


and  therefore,  using  these  bounds  in  (6.2),  we  obtain 

2hJ  max  j-  |J?'(x)| 

6  - 1|  < - rr - • 
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Because  of  the  assumed  smoothness  of  u  and  the  regularity  of  the  mesh  (i.e.  0  <  C\  < 
X^7  <  Cj  <  oo),  /P(x)  =  C?(hf),  for  x  €  therefore,  whenever  |e.|,  ^  >  Ch^,  (C  >  0), 
then  |6  —  1|  — *  0  as  h,  -*  0;  that  is,  B-M  indicator  is  asymptotically  exact.  Q 

Let  us  underline  the  role  of  the  smcthness  of  u.  In  the  proof  of  the  theorem  we  assume 
that  u"'(x)  is  bounded  for  x  €  Instead,  it  could  be  assumed  u  €  H2+<(I,),  («  >  0). 

2h$  max  ~  |R'(t)| 

The  elemental  effectivity  index  depends  on  the  quotient  Q,  := - [7717^ - and  the 

asymptotic  exactness  occurs  if  Qi  — *  0  as  h,  — *  0.  Nevertheless,  for  fixed  ht,  depends 
on  Q ,  and  a  high  quality  index  requires  Qi  <  1.  We  shall  return  to  this  problem  in  the 
subsequent  sections. 
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THEOREM  6.2.  Z-2  indicator  is  not  asymptotically  exact. 

Proof.  Assume  that  u  is  a  quadratic  polynomial;  consider  a  mesh  such  that  x  = 
hi+ 1  =  2 h{.  An  explicit  computation  shows  that  the  elemental  effectivity  index  is  f,  =  2 
independently  of  the  meshsize  ht.  Therefore  Z-Z  is  not  asymptotically  exact.  0 

Let  us  now  consider  as  a  numerical  example  a  particular  problem  (4.1)  whose  solution 
is  given  by  u(x)  =s  9l(xp+t^).  Z-Z  and  B-M  indicators  are  identical  for  uniform  meshes, 
hence  to  observe  different  performances  we  need  to  use  a  non  uniform  one.  Let  x,_2  =  0.01, 
Xi_i  =  0.04,  Xi  =  0.09,  Xi+i  =  0.16.  Table  6.1  shows  the  elemental  effectivity  index  if,  for 
both  indicators.  _ 


?  = 

0.0 

q  = 

0.5 

p 

Z-Z 

B-M 

Z-Z 

B-M 

0.55 

1.07 

1.78 

4.70 

6.19 

0.65 

1.03 

1.68 

3.18 

5.21 

0.75 

0.99 

1.57 

1.59 

2.80 

0.85 

0.97 

1.48 

1.18 

2.07 

0.95 

0.96 

1.39 

1.01 

1.69 

1.05 

0.96 

1.32 

0.94 

1.45 

1.15 

0.96 

1.26 

0.92 

1.29 

1.25 

0.97 

1.20 

0.93 

1.17 

1.35 

0.98 

1.15 

0.97 

1.08 

1.45 

1.01 

1.11 

1.01 

1.02 

1.55 

1.03 

1.07 

1.07 

0.99 

1.65 

1.07 

1.04 

1.13 

0.97 

1.75 

1.10 

1.02 

1.20 

0.96 

1.85 

1.15 

1.01 

1.27 

0.97 

1.95 

1.19 

1.00 

1.35 

0.99 

2.05 

1.24 

1.00 

1.44 

1.02 

2.15 

1.30 

1.01 

1.52 

1.06 

2.25 

1.35 

1.02 

1.61 

1.11 

2.35 

1.41 

1.04 

1.71 

1.16 

2.45 

1.48 

1.06 

1  81 

1.21 

Table  6.1.  Elemental  effectivity  indices  for  different  singular  solutions. 

We  observe  that  the  Z-Z  indicator  is  better  for  nonsmooth  solutions  ( p  <  1.5)  and  the 
B-M  indicator  is  better  for  smooth  u.  The  inaccuracy  of  the  B-M  indicator  for  p  small  is 
directly  related  to  the  nonsmoothness  of  the  function  u  and  the  poor  quality  of  the  Z-Z 
indicator  for  p  large  is  related  to  the  fact  that  it  is  not  asymptotically  exact.  Table  6.1 
shows  that  to  measure  the  quality  of  an  indicator  is  not  a  simple  task.  In  the  next  sections 
we  shall  address  in  more  detail  the  observed  effects. 

Let  us  consider  now  a  particular  problem  (4.1)  whose  solution  is  given  by 
(6.3)  u(x)  =  i*  —  x,  x  €  /  =  (0, 1)  . 
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For  p  >  1.5,  u  €  H2+t(I),  (<  >  0),  and  hence  the  asymptotic  exactness  of  the  B-M  indicator 
should  be  reflected  in  its  performance.  But  for  p  <  1.5  we  cannot  expect  to  observe  high 
quality  effectivity  indices  for  this  indicator. 

Remark  6.1.  We  study  these  kind  of  problems  with  solutions  (6.3)  since  they  are  one 
dimensional  models  for  the  singularity  of  the  solutions  of  two  dimensional  elliptic  problems 
in  the  neighborhood  of  a  comer.  □ 


We  consider  graded  meshes  with  nodes 


(6.4) 


=  0,...,n 


with  0  >  0.  Table  6.2  shows  the  true  error  and  the  elemental  effectivity  indices  for 
n  =  10,  for  p  bb  1.25  and  p  =  2.25,  and  for  0  =  2.0  and  0  —  0.5. 


P  = 

1.25 

II 

to 

o 

0  SB  0.5 

Int. 

error 

Z-Z 

B-M 

error 

Z-Z 

B-M 

1 

6.45(-3) 

1.172 

0.391 

8.61(-2) 

0.256 

0.618 

2 

7.93(-3) 

0.937 

1.441 

8.85(-3) 

3.574 

1.392 

3 

8.02(-3) 

0.968 

1.198 

4.86(-3) 

1.254 

1.034 

4 

8.04(-3) 

0.983 

1.103 

3.31(-3) 

1.116 

1.016 

5 

8.05(-3) 

0.990 

1.063 

2.49(-3) 

1.068 

1.009 

6 

8.06(-3) 

0.993 

1.042 

1.98(-3) 

1.044 

1.006 

7 

8.06(-3) 

0.995 

1.030 

1.64(-3) 

1.032 

1.004 

8 

8.06(-3) 

0.996 

1.023 

1.40(-3) 

1.024 

1.003 

9 

8.06(-3) 

0.997 

1.018 

1.21(-3) 

1.018 

1.002 

10 

8.07(-3) 

0.970 

1.084 

1.07(-3) 

1.080 

1.021 

P  = 

2.25 

0  —  2.0 

0  =  0.5 

Int. 

error 

Z-Z 

B-M 

error 

Z-Z 

B-M 

1 

2.11(-4) 

4.024 

1.341 

8.91(-2) 

0.469 

1.131 

2 

1.67(-3) 

1.807 

1.050 

3.02(-2) 

1.879 

0.969 

3 

4.57(-3) 

1.352 

1.019 

2.17(-2) 

1.116 

0.996 

4 

8.93(-3) 

1.194 

1.010 

1.75(-2) 

1.054 

0.998 

5 

1.47(-2) 

1.122 

1.006 

1.50(-2) 

1.031 

0.999 

6 

2.20(2) 

1.083 

1.004 

1.32(-2) 

1.021 

0.999 

7 

3.07(-2) 

1.060 

1.003 

1.19(-2) 

1.015 

0.999 

8 

4.09(-2) 

1.045 

1.002 

1.09(-2) 

1.011 

1.000 

9 

5.25(-2) 

1.036 

1.002 

1.01(-2) 

1.008 

1.000 

10 

6.56(-2) 

0.872 

0.974 

9.38(-3) 

1.050 

0.993 

Table  6.2.  Elemental  effectivity  indices  for  different  graded  meshes. 


20 


We  observe  the  following  features: 

i)  Z-Z  indicator  is  of  high  quality  when  the  mesh  is  equilibrated;  i.e.  the  mesh  is 
such  that  the  errors  are  almost  equidistributed  in  all  the  subintervals.  (In  the  next 
section  we  shall  address  this  feature  in  detail). 

ii)  B-M  indicator  is  in  general  of  higher  quality,  except  for  those  meshes  that  are  well 
equilibrated  and  when  the  solution  is  not  sufficiently  smooth. 

iii)  The  performance  of  Z-Z  and  B-M  indicators  in  the  first  and  the  last  intervals  is  in 
general  poor.  This  feature  is  related  to  extension  aspects. 

7.  Further  analysis  of  Z-Z  and  B-M  indicators  and  their  generalization.  As 
we  have  seen  in  section  4,  U(x, )  is  a  weighted  average  of  u'pE  on  I,  and  7,+  j .  For  ooth 
indicators,  Z-Z  and  B-M,  U(x,)  lies  between  u'n | /,  and  u'pJ/i+l ;  in  fact,  according  to  (4.4), 
(4.11)  and  (4.12): 


^(*0  =  u'f,Ia  +  <*■  KJ,  *  “«Ia  +  <*•  (“feU+i  -««!/.■)  - 


with 

(7.1) 

for  Z-Z  indicator  and 

(7.2) 


1 

1  +  A,  ’ 


a. 


a 


BM 


1  +  Afl  ■ 


for  B-M  indicator,  where  A,  :* 

Let  us  now  address  the  problem  of  how  to  select  a*.  Using  (4.7),  (4.8)  and  the  relation 
a,  =  1  +  0i,  the  indicator  can  be  written: 


(7.3)  7?  -  y  [(l-a,_i)J  K,]*,!  -°i(l  -a«-, )Ke]-i  KJ,  +  q?  Ktl! 


For  x  €  Ii  let 


(7.4)  u'(x)  =  u'p,!a  +  al(x-xi_1/2)  +  r,(x)  , 

with  x,_i/j  :=  2izX±2i.  and 

o' yj*  -  =  e-(l!(I  -  Xi-'n)dx  ~  -if/,,  e'(x)dx  ■ 

hence  r\(x)  is  orthogonal  to  1  and  to  (x  —  x;_i^),  (that  is,  Jr  r,(x)dx  =  0  and  r,(x)(x- 

T>-\fl)dx  =  0). 
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Let  us  assume  that  u"(x)  >  0  (or  u"(x)  <  0)  for  all  i  6  J*.  Then  e«(x)  <  0  (resp. 
e«(x)  >  0)  Vx  €  (xj-i,Xj),j  =  »-l,i,i  +  l,  and  *o  a<_i , a,, at+1  >  0  (resp.  ai_i,a,,a,+1  < 
0).  Let  L,  and  R,  be  defined  by 


it/  \  ai^i  .  /  ajh>  /,  T  \ 

= - 2~  +  - - 2~^  ~  Li ^ 


(7.5) 


=  ^  +  r,(xt)  =  ^(l  +  *,)• 


Let  K,  be  defined  by 


(7-6)  |«.| 


2  f  2/  -j  t'2\ 


If  u  is  sufficiently  smooth,  Li,Rt,Kt  =  0(/ij). 

Under  the  assumption  u"(x)  >  0  (or  u"(x)  <  0)  Vx  €  /<,  the  jumps  and 

[u'rE]  are  strictly  positive  (resp.  strictly  negative)  and  hence  there  exist  coefficients  o,_j 
and  ai  such  that 


(7.7) 


Using  these  coefficients  in  (7.3)  we  obtain  for  the  corresponding  estimator 

that  is,  the  coefficients  defined  in  (7.7)  are  optimal  in  the  sense  that  they  yield  an  exact 
estimator  of  the  error. 

By  using  (7.5)  the  optimal  coefficient  a i  can  be  written 


(7.8) 


a.  = 


(W) 


1  +  (ttT)  i1^) 

and  using  (7.6),  since  a,  and  a,+  i  have  the  same  sign, 


(7.9) 


a,  = 


(^) 

1  +  1 

(ttT) 

l( 

.  »+*•?  ) 
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A  similar  expression  holds  for  the  other  optimal  coefficient  a,_i . 

Therefore,  whenever  the  mesh  is  equilibrated  with  respect  to  the  energy  norm  (i.e. 
85  25  le»li,/,+1 )»  u  13  smooth  enough  as  L„  Rt,  L,+i ,  K,  =  o(  1)  to 

hold,  then  the  optimal  coefficients  are 


a.-i  = 


1 


1  +  A 


* 


+  o(l) 


and 


a,  - 


i-i 


1  +  A* 


+  o(l) 


and  we  have  the  following  theorem. 


THEOREM  7.1.  Let  u"(x)  >  0  (or  u"(x)  <  0 )  for  edl  x  T  and  let  R,,L,,K,  =  o(l), 
t  =  l,...,n.  Then,  if  the  mesh  is  equilibrated  with  respect  to  the  energy  norm,  the 
coefficients 


(7.10) 


1 

1  + A? 


yield  an  asymptotically  exact  estimator. 


The  values  of  L,,  Rt  and  K,  govern  the  effectivity  of  the  indicators.  If  the  solution  u 
is  less  smooth,  the  values  of  these  constants  are  larger  and  we  have  to  expect  less  accurate 
effectivity  indices.  Wj  shall  show  that  when  the  solution  has  a  very  strong  singularity,  the 
estimator  defined  by  a°pt  is  not  necessarily  optimal. 

Consider  the  model  probiem  of  section  6  with  solution  (6.3).  In  this  case  u"( r)  does 
not  change  sign  in  I  and  so  Theorem  7.1  can  be  applied.  However  for  £  <  p  <  1,  Li,  R, 
and  K,  are  not  negligible  and  become  larger  as  p  — ►  In  fact,  in  this  particular  case,  when 
using  meshes  equilibrated  with  respect  to  the  energy  norm,  we  shall  see  that  choosing 


1 

1  +  A? 


with  7  >  £  can  be  preferable. 

In  Table  7.1  we  show  the  result'  for  various  values  of  p  and  the  corresponding  values 
of  j3  in  (6.4)  leading  to  asymptotically  equilibrated  meshes. 
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0.75  -  0 


7  =  1 


1.05564 

1.00429 

0.99922 

0.98554 

0.98309 

0.98470 

0.98708 

0.98925 

0.99106 

0.83594 


0.95241 

0.75513 

0.86071 

0.92311 

0.95261 

0.96809 

0.97710 

0.98278 

0.98659 

0.98765 


=  1.25  -  0  =  2 


7  =  1 

7  =  \ 

1.17157 

0.99033 

0.93675 

0.90257 

0.96772 

0.97650 

0.98300 

0.98921 

0.98962 

0.99374 

0.99302 

0.99590 

0.99499 

0.99710 

0.99623 

0.997S4 

0.99707 

0.99833 

0.96997 

0.99S4S 

=  1.75  -  0  = 


error 

7  =  1 

7  =  1 

1.05414 

0.99399 

jUmi  7KXW 

0.95943 

0.97301 

m/m  fKSI 

0.99164 

0.99574 

1.57(-2) 

0.99608 

0.99809 

1.57(-2) 

0.99770 

0.99890 

1.57(-2) 

0.99848 

0.99928 

1.57(-2) 

0.99892 

0.99949 

1.57(-2) 

0.99920 

0.99962 

1.57(-2) 

0.99938 

0.99971 

1.57(-2) 

0.99414 

0.99974 

=  2.25  -  0 


7  =  1 


0.95182 

1.04344 

1.00638 

1.00292 

1.00169 

1.00111 

1.00079 

1.00059 

1.00046 

1.00414 


1.00696 

1.02352 

1.002SS 

1.00128 

1.00074 

1.0004S 

1.00034 

1.00025 

1.00020 

1.0001S 


Table  7.1.  Comparison  of  indicators  for  equilibrated  meshes  with  respect  to  the  energy  norm. 

This  table  shows  that  for  p  =  .75,  7  =  1  (i.e.  Z-Z  indicator)  yields  better  effectivity 
indices  than  7  =  i.  The  particular  good  performance  of  this  estimator  is  related  with 
the  monotonicity  of  u  and  is  not  always  valid  when  the  solution  is  highly  singular  (see 
Table  6.1).  For  p  =  0.75  and  for  p  —  1.25,  the  effectivity  index  in  the  first  elements  of  the 
estimator  corresponding  to  7  =  j  is  not  so  close  to  one,  because  the  mesh  is  not  so  well 
equilibrated. 

Remark  7.1.  Assume  once  more  that  u  is  sufficiently  smooth  for  L ,,  R,  and  I\,  to 
be  negligible  in  (7.5)  and  (7.6).  Then  the  optimal  coefficients  (7.8)  are 

(7.11)  a*  - - t -  . 

1  + 

<*i  fit 

!•«  Kill.,;,)  =”  [/,,  |e,,(x)|*<fij  *  be  the  usual  L3(I,)~ norm  of  e'u.  A  simple  Taylor  series 
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argument  applied  to  (7.4)  shows  that  if  u  is  smooth  enough  then 


If  the  mesh  is  equilibrated  with  respect  to  this  norm,  then  the  optimal  coefficients  are 


(7-12) 


an  ss  a(1*)  := 


1  + A* 


up  to  higher  order  terms.  Obviously  (7.10)  is  a  particular  case  of  this  formula  with  s  =  2. 
But  for  instance,  if  the  mesh  is  equilibrated  with  respect  to  the  norm  ||e„ II £, i ( /  )>  then  Z-Z 
indicator  is  the  optimal.  On  the  other  hand,  when  s  — »  oo,  a*,)  — ♦  i;  that  is,  the  simple 
average  of  u^l/,  and  u'FE|/,+l  is  the  correct  value  of  U(x,)  when  the  mesh  is  equilibrated 
with  respect  to  the  norm  ||e'J|  £,<»(/;)•  D 

If  we  knew  the  exact  errors  jeu  | }  j  for  a  given  mesh,  then  we  could  use  them  in  (7.9) 
to  obtain  optimal  coefficients 


(7.13) 


Using  this  information,  we  could  improve  the  effectivity  indices.  We  show  this  effect  in  the 
following  table. 


1 

P 

=  0.75  - 

3  =  9.0 

P 

=  1.75  - 

3  =  2.5 

Int. 

error 

o'," 

<•!" 

q(2) 

« 

error 

,Ji) 

a. 

a!2> 

s;!i 

B 

BEB 

1.095 

■SSI 

3.247 

■SOBS 

1.118 

mm 

1.011 

Ifffl 

1.656 

1 

BBBM 

1.128 

mSm 

0.907 

WEfim 

1.2S6 

1.146 

1.95C-2) 

1.110 

in 

0.929 

8.59(-3) 

1.156 

1.076 

Epilfl 

2.42(-2) 

1.084 

mm 

0.950 

1.29(-2) 

1.097 

1.047 

2.86(-2) 

1.063 

mm 

0.964 

1.79(-2) 

1.066 

1.031 

3.27(-2) 

0.961 

0.973 

2.35(-2) 

1  04S 

1.023 

3.66(-2) 

0.970 

0.979 

2.97(-2) 

1.036 

1.017 

1.000 

4.04(-2) 

0.976 

0.984 

3.64(-2) 

1.02S 

1.013 

1.000 

10 

4.40(-2) 

0.700 

0.936 

0.986 

4.36(-2) 

0.S75 

0.914 

1.000 

Table  7.2 
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Obviously  the  exact  errors  |e„|j  r  are  not  available.  Nevertheless,  if  we  use  the  esti¬ 
mated  errors  and  if  their  ratio  is  reasonably  close  to  the  ratio  of  the  true  errors,  then  it  is 
possible  to  iterate.  Table  7.3  shows  the  results  obtained  by  this  iterative  procedure. 


P 

=  0.75  - 

(3  =  9.0 

P 

=  1.75  - 

(3  =  2.5 

Int. 

error 

o<2) 

Iter. 

s<,2) 

error 

«i2) 

Iter. 

~(2) 

1.99(-3) 

1.237 

wgm 

1.095 

3.56(-4) 

ppa 

8. 1 1(-3) 

0.881 

ESI 

1.011 

2.15(-3) 

cm 

'  tM 

1.42(-2) 

0.S61 

HI 

0.907 

4.96(-3) 

1.146 

1.034 

1.95(-2) 

0.897 

0.934 

0.929 

8.59(-3) 

1.076 

0.980 

l  | 

2.42(-2) 

0.928 

0.957 

0.950 

1.29(-2) 

mm 

1.019 

1.001 

2.86(-2) 

0.948 

0.962 

0.964 

1.79(-2) 

mm 

0.986 

1.001 

7 

3.27(-2) 

0.961 

0.981 

0.973 

2.35(-2) 

BUI 

1.013 

1.001 

8 

3.66(-2) 

0.970 

0.971 

0.979 

2.97(-2) 

1.017 

0.989 

1.000 

9 

SEE! 

0.976 

0.998 

0.984 

3.64(-2) 

1.013 

1.010 

1.000 

10 

dsn 

0.936 

0.968 

0.986 

4.36(-2) 

0.914 

0.991 

1.00G 

Table  7.3 

As  a  conclusion,  let  us  remark  that  for  an  adaptive  procedure  based  on  equilibrated 
meshes,  it  is  convenient  to  use  the  proper  optimal  indicator.  The  Z-Z  indicator  is  not 
the  proper  one  when  the  adaptive  procedure  is  based  on  equilibration  with  respect  to  the 
energy  norm.  Nevertheless  it  is  reasonably  close,  but  the  proper  one  is  preferable.  The 
estimates  could  be  improved  if  a-posteriori  information  about  the  errors  is  used,  namely 
the  quotient  of  computed  estimated  errors.  The  results  for  the  first  and  last  intervals  are 
less  reliable  because  of  extension  aspects. 

8.  Robustness  and  quality  measure  of  an  indicator.  In  section  4  we  introduced 
a  family  of  indicators  based  on  coefficients  a,,  j  (see  (4.7-8))  and  in  the  subsequent 
sections  we  discussed  their  performance  from  an  asymptotic  point  of  view.  This  approach 
cannot  be  well  applied  when  the  asymptotic  assumptions  (i.e.  Li,  R,,  I\,  =  o{l)  as  in 
section  7)  are  not  satisfied. 

In  (5.4)  we  introduced  another  measure  of  the  quality  of  an  indicator: 


||£,||  :=  sup 


K  -  VM  o./. 


that  can  be  used  when  the  solution  is  known  to  belong  to  a  subspace  H  C  Hq(I).  Now, 
based  on  this  measure,  we  shall  develop  principles  for  judging  an  indicator  without  asymp¬ 
totic  assumptions.  (From  now  on,  we  write  V,u,  V,v  instead  of  V,  in  (4.7),  since  we  shall 
need  to  deal  with  functions  V,  defined  by  solutions  u,  u  of  different  problems). 
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Since  the  indicators  defined  by  (4.7-8)  satisfy  |e'„  —  Viu|0  j.  =*  0  when  v  is  a  linear 
function  in  then  we  can  analyze  the  space  'H  modulo  V\  (the  space  of  linear  functions). 
One  of  the  possibilities  is  to  consider  H2(Ii)/V i;  however,  let  us  note  that  ||£,|j  is  well 
defined  even  for  Hl(Ii)/V\. 

On  the  other  hand,  £,  depends  on  the  particular  indicator  (i.e.  on  a,  and  /3,_i  in  (4.7)), 
on  the  mesh  T  and  on  the  particular  space  Ti..  We  shall  write  explicitly  this  dependence; 
so,  instead  of  (5.4),  let  us  call 


(8.1) 


7 l(tt,a,0,7) 


sup 


K  -  'Mo./,. 
HI* 


sup 

v€ft/?i 


K-vMom 

IMIft/-?! 


the  robustness  index  of  the  indicator;  (from  now  on  we  omit  the  subindices  of  the  coeffi¬ 
cients  a,  and  /3,_i  defining  Vju). 

If  we  do  not  know  anything  more  about  the  solution  u  than  that  u  £  H,  then  the 
safest  choice  of  the  indicator  is  one  whose  robustness  index  is  minimal.  In  fact,  for  any 
pair  (a,  (3)  there  exists  u  €  H  such  that 


7 Z(H,a,0,T) 


K  -  vtf|0>/i 


n?< 


and  so  we  accept  the  possibility  that  u  can  be  the  solution  of  our  problem.  Nevertheless,  u 
(which  depends  on  H)  could  be  an  unlikely  candidate  as  the  solution  of  the  problem  under 
consideration.  This  would  indicate  that  the  space  H  is  too  large  for  our  purpose.  So,  the 
role  of  7 i  for  the  quality  assesment  is  obvious. 

Sometimes,  it  is  known  a-priori  that  the  solution  u  belongs  to  all  the  spaces  of  a 
certain  family.  To  deal  with  this  case,  we  introduce  some  notation.  Let  ACr,*  denote  a 
one-parameter  family  of  spaces  7f7,  7  €  T,  where  i/>  :  T  — ►  R  is  a  positive  function;  we  say 
that  u  €  ACr,*>  if  u  €  H7,  V7  €  T  and 


(8.2)  3C  >  0  :  0(7)IM|*,  <  C  ,  V7  €  T  . 

For  example,  if  T  =  N,  Hn  =  Hn(I)  and  i/>(n)  =  n!e“n,  n  €  N,  then  (8.2)  holds  if  and  only 
if  u  is  analytic  in  /;  in  [2,3] ,  these  kind  of  spaces  have  been  used  for  the  characterization  of 
the  regularity  of  the  solutions  of  elliptic  partial  differential  equations  in  the  neighborhood 
of  a  corner  of  the  domain. 

If  for  a  given  problem,  we  know  not  only  that  ti  €  ACr.*-,  but  also  Hull*,/?*,  for  7  €  T, 
we  could  then  compute 

(8.3)  G(ACr.*,«,/J,T,u):«  inf  {7J(W7,a,/?,r)||u||ww)>|  } 

and  since  \e'v  -  Vju|0  It  <  7Z(7i^,  a,  7’)!l«!|*.i/-pl  for  all  7  6  T,  Q{Kr^,a,/3,T,u)  could 
then  be  us<*d  to  asses  the  quality  of  the  indicator  for  this  particular  problem. 
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In  practice,  we  may  know  only  approximately.  If  we  denote  the 

approximation  of  \\u\\-n^ /-px  obtained  by  using  the  (a-posteriori)  information  that  we  have 
about  the  solution,  then  we  call 

(8.4)  <2*(£l\*,  a,  0,  r,  u)  :=  inf  { a,  /?,  T )||ti||^ /t>i  } 

the  quality  index  of  the  indicator  for  this  particular  problem. 

9.  Computation  of  robustness  and  quality  indices.  In  this  section  we  shall 
address  as  an  example  the  computation  of  the  indices  defined  above  for  some  weighted 
Sobolev  spaces  modeling  the  space  of  solutions  of  elliptic  PDEs  in  the  neigboi-kood  of 
comers. 

Let  /  be  an  interval  and  m  >  /  >  0  be  integers,  then  by  we  denote  the 

weighted  Sobolev  space  that  is  the  completion  of  the  set  of  the  infinitely  differentiable 
functions  under  the  norm 


:= 

7 ,1 


l-l 


E  il»“>(*)il,.r+  E  ll*’"*''  »“>(*)llo,r 

.icO  *~J 

where  *s  vo‘d  and  :=  For  Jfc  =  »,...,  m,  let 

\v\(y  :=  u(*)(i)||0j- 


1/2 


Let  us  introduce  first  some  auxiliary  theoretical  results. 

LEMMA  9.1.  If0<7<7'<l  and  m  >  /  >  0,  then  the  inclusion  ff™+l,l(f)  «— » 
Hy'l(I)  is  compact. 

Proof.  See  (21)  (page  287).  Q 

Theorem  9.1.  Let  (a,/?)  e  R2  and  7  €  (0,1),  then 


(9-1) 


K  -  V,v\ 


i.e.  (e'u  —  V,v)  is  a  bounded  linear  operator  of  H\2'2\l,)/V  1  into  L2(I,). 

Proof  Since  the  inclusion  H2,2(I )  •— »  Hl(I)  is  continuous,  then  (9.1)  holds  when  on 
the  right  hand  3ide  is  used  instead  of  llullH>.»(J')/p1  •  Since  e'.  -  V,u  =  0  for 

v  6  Pi,  (9.1)  follows  immediately.  Q 

So  we  have: 
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COROLLARY  9.1.  The  robustness  index  is  defined  for  any  (a,/3)  €  R2,  any  mesh  T 
and  H  —  H2'2(Ii),  with  0  <  7  <  1. 

For  0  <  7  <  1,  let 

(v,  u?)2'2  :*  f  [x2’,+2t/"(x)u>"'(x)  +  x27v"(x)u;"(x)]  dx 

JZ 

denote  the  scalar  product  in  the  Hilbert  space  H2'2{Ii)/V\  and  |  •  |2,2  :=  (•,-)2’2  ^  the 
corresponding  norm. 

For  any  v  €  H2,2(Ii)/V\,  there  exists  a  unique  zv  €  H2'2(It)/V  1  such  that 
(9.2)  =  f[  («'„  -  Vtv)(e'w  -  V,w)dx  ,  Vu,  €  H2'2(I,)/P1  . 

(this  follows  from  Theorem  9.1).  Hence,  T(v)  :=s  zv  defines  an  operator  of  H2-2(I,)/Vi 
into  itself  and  we  have  the  following  theorem. 

THEOREM  9.2.  T  is  a  self-adjoint,  positive  definite,  compact,  linear  operator  of  H2,2 (I X)/V\ 
into  itself. 

Proof.  It  is  a  simple  consequence  of  the  definition  and  Lemma  9.1.  Q 
Thus  we  have: 

Corollary  9.2. 

K(H2'2(T,),a,0,T)  =  A*  , 

where  X  is  the  maximal  eigenvalue  of  the  problem  T(  v)  =  Xv. 

<2 

Let  us  now  consider  the  mesh  T  with  nodes  x,  :=  (^)  ,  (/,  :=  (x,_i ,  x,)).  In  Table  9.1 
we  show  TZ(H2,2(It),  a,  0,  T)  for  7  =  0.0  and  7  =  0.9  at  each  interior  subinterval,  for  Z-Z 
and  B-M  indicators  and  for  the  optimal  indicator  in  the  sense  of  minimising  V..  7  =  0.0 
corresponds  to  a  space  of  smooth  functions  and  7  =  0.9  to  a  set  of  functions  with  a  strong 
singularity  at  the  origin. 

We  observe  that  for  very  nonsmooth  functions  (7  =  0.9)  the  robustness  index  indicates 
that  the  Z-Z  indicator  performs  better  than  the  B-M  mdicator  on  those  intervals  closest 
to  the  singularity  (i  =  2,3,4).  On  the  other  hand,  for  smooth  functions  (7  =  0.0)  B-M 
robustness  indices  are  always  better. 

Now  let  £r>  be  the  one-parameter  family  of  spaces  H2,2(Ij),  7  6  T  :=  (0, 1),  with 
0(7)  s=  1,  V7  €  F,  (I2  =  (0.01,0.04)  and  I?  =  (0.00,0.09)).  In  Table  9.2  we  show  the  quality 
indices  Q(ACr,*,£>,/3,T,u)  for  problems  which  solutions  are  the  functions  u(x)  =  ■R(x,’+,, ) 
for  various  values  of  p  and  q  =  0  0  and  q  =  0.5.  We  use  the  exact  value  of  Hull  5,r  , to 
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compute  Q.  We  also  show  in  this  table,  the  exact  error  |e'B|0  /j,  the  estimated  error  |Vju|0/i, 
the  effectivity  index  £ 2 >  the  error  in  the  estimate  of  e'„  (i.e.  err.  est.:=  |e'„  —  V,u|0  and 

\t  —  V,  u  j 

the  percentage  of  this  last  error  (i.e.  %:=  100  *|.n -■* ). 


7  =  0.0 

7  =  0.9 

Int. 

Z-Z 

B-M 

optimal 

Z-Z 

B-M 

optimal 

0.451(-4) 

0.232(-4) 

0.103(-4) 

0.139(-1) 

0.645(-l) 

0.897(-2) 

0.630(-4) 

0.363(-4) 

0.165(-4) 

0.542(-2) 

0.865(-2) 

0.254(-2) 

0.731(-4) 

0.440(-4) 

0.229(-4) 

0.227(-2) 

0.244(-2) 

0.983(-3) 

5 

0.786(-4) 

0.484(-4) 

0.274(-4) 

0. 1 10(-2) 

0.978(-3) 

0.481(-3) 

6 

0.819(-4) 

0.509(-4) 

0.307(-4) 

0.595(-3) 

0.476(-3) 

0.269(-3) 

7 

0.839(-4) 

0.525(-4) 

0.336(-4) 

0.348(-3) 

0.262(-3) 

0.160(-3) 

8 

0.852(-4) 

0.536(-4) 

0.363(-4) 

0.218(-3) 

0.157(-3) 

0.l02(-3) 

9 

0.861(-4) 

0.543(-4) 

0.378(-4) 

0.143(-3) 

0.101(-3) 

0.685(-4) 

Table  9.1.  Robustness  indices. 


<1 

=  0.0 

Z-Z 

P 

Klo./, 

f • 

0" 

6 

err.  est. 

% 

<2 

0.55 

0.94(-l) 

0.11(-0) 

1.14 

0.52(-l) 

55.6 

0.21(-0) 

0.75 

0.33(-l) 

0.29(-l) 

0.84(-2) 

25.8 

0.40(-1) 

0.95 

0.32(-2) 

0.11(-2) 

28.6 

0.42(-2) 

1.15 

mESSm 

0.57(-2) 

0.89 

0.28(-2) 

43.1 

0.66(-2) 

1.35 

0.82(-2) 

0.81(-2) 

0.99 

0.48(-2) 

58.7 

0.84(-2) 

1.55 

0.69(-2) 

1 mmm 

1.13 

0.52(-2) 

74.8 

0.73(-2) 

1.75 

0.50(-2) 

1.29 

0.46(-2) 

91. S 

0.56(-2) 

1.95 

0.33(-2) 

0.49(-2) 

1.48 

0.37(-2) 

110.5 

0.42(-2) 

2.15 

0.21(-2) 

mm 1 

1.69 

0.28(-2) 

131.3 

0.31(-2) 

2.35 

0. 13( -2) 

1.93 

154.8 

0.22(-2) 

Table  9.2.a.  Quality  indices  Q  of  Z-Z  indicator  for  u(x)  =  3?(xp+,,),  ( q  =  0.0). 
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Table  9.2.d,  Quality  indices  Q  of  B-M  indicator  for  u(i)  =  ( q  =  0.5). 

Since,  in  practice,  the  exact  norm  of  the  solution  is  not  known,  we  need  to  estimate  it. 
The  second  derivative  can  be  recovered  from  the  error  estimate  and  the  third  one  by  using 
differences  of  these  estimates  on  neighboring  elements.  By  this  way  we  can  compute  the 
quality  indices  <2*.  In  Table  9.3  we  report  the  values  of  Q"  corresponding  to  tables  9.2. 


9  = 

0.0 

q  = 

0.5 

? 

Z-Z 

B-M 

Z-Z 

B-M 

0.55 

0.12(-0) 

0.19(-0) 

0.13(-0) 

0.25(-0) 

0.75 

0.38(-l) 

0.51(-1) 

0.76(-l) 

0.11(-0) 

0.95 

0.42(-2) 

0.49(-2) 

0.43(-l) 

0.47(-l) 

1.15 

0.68(-2) 

0.71(-2) 

0.23(-l) 

0.21(-1) 

1.35 

0.87(-2) 

0.81(-2) 

0.13(-1) 

0.99(-2) 

1.55 

0.75(-2) 

0.64(-2) 

0.76(-2) 

0.55(-2) 

1.75 

0.58(-2) 

0.45(-2) 

0.51(-2) 

0.37(-2) 

1.95 

0.43(-2) 

0.32(-2) 

0.41(-2) 

0.29(-2) 

2.15 

0.31(-2) 

0.22(-2) 

0.37(-2) 

0.27(-2) 

2.35 

0.23(-2) 

0.17(-2) 

0.36(-2) 

0.26(-2) 

Table  9.3.  Quality  indices  Q*. 

From  tables  9.2  and  9.3  we  observe  the  following  features. 

1)  Sometimes,  although  V,u  does  not  approximate  the  error  e'B  well,  ||  V, u ||0  ^  can  be 
a  close  approximation  of  ||e'J|0 

2)  Z-Z  indicator  performs  better  than  B-M  when  the  solution  is  nonsmooth  (provid¬ 
ed  that  the  mesh  is  properly  graded)  and  B-M  performs  better  if  the  solution  is 
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smooth. 

3)  The  quality  indices  correctly  indicate  that  Z-Z  indicator  is  preferable  for  nonsmooth 
solutions  while  B-M  is  preferable  for  smooth  ones.  In  principle,  a  criterion  for 
selecting  Z-Z,  B-M  or  any  other  indicator  could  be  based  on  comparing  their  quality 
indices. 

10.  Conclusions.  In  the  one-dimensional  setting  we  have  shown: 

a)  From  an  asymptotic  point  of  view,  an  optimal  indicator  can  be  derived  for  any 
adaptively  constructed  mesh.  The  Z-Z  indicator  is  optimal  when  the  adaptive 
procedure  equilibrates  the  elemental  error  in  the  I.1 -norm  and  is  not  optimal  if 
the  energy  norm  (i.e.  £2-norm  of  e',)  is  equilibrated. 

b)  Z-Z  indicator  is  not  asymptotically  exact  for  a  general  mesh  while  B-M  is. 

c)  In  a  nonasymptotic  sense,  when  the  solution  is  nonsmooth  and  the  mesh  is  reason¬ 
ably  graded,  Z-Z  indicator  gives  better  results  than  B-M.  Otherwise  B-M  does. 

d)  The  quality  and  robustness  indices  can  be  used  to  select  am  indicator. 

In  the  two-  dimensional  setting  we  have  shown: 

a)  For  general  meshes,  an  asymptoticadly  exact  estimator  cam  not  be  achieved  by  a 
simple  L2  -projection  technique. 

b)  For  translation  invariant  meshes,  asymptotically  exact  estimators  cam  be  defined 
by  utilizing  superconvergence  techniques  based  on  differences. 

c)  For  more  general  meshes,  the  way  to  define  a  reasonable  estimator  is  very  likely  to 
assume  that  the  mesh  is  locally  translation  invariamt  (ad though  it  is  not). 

d)  The  major  principles  analyzed  in  the  one-dimensional  setting  could  be  utilized  in 
two  dimensions,  although  further  research  work  is  needed  in  this  direction. 
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those  in  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  includes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology. 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agencies 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  BabuSka.  Chairman.  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Science  and  Technology,  University  of  Maryland,  Colleee 
Park,  Maryland  20742-2431. 


